Mimicking a Squeezed Bath Interaction: 
Quantum Reservoir Engineering with Àtoms 



N. Lütkenhaus 
J.I. Cirac 
P. Zoller 

■ Institut für Theoretische Physik 



O 

o 

On 
Oi. 



Universitat Innsbruck 



Abstract 



ç-C) . The interaction of an atòmic two-level system and a squeezed vacuum leads 

o, 

to interesting novel effects in atòmic dynamics, including line narrowing in 



resonance fluorescence and absorption spectra, and a suppressed (enhanced) 
decay of the in-phase and out-of phase component of the atòmic polarization. 
On the experimental side these predictions have so far eluded observation, 



essentially due to the difficulty of embedding àtoms in a 4tt squeezed vacuum. 
In this paper we show how to "engineer" a squeezed-bath-type interaction for 
an effective two-level system. In the simplest example, our two-level atom 
is represented by the two ground levels of an atom with angular momentum 
J = l/2— >J = l/2 transition (a four level system) which is driven by 
(weak) làser fields and coupled to the vacuum reservoir of radiation modes. 
Interference between the spontaneous emission channels in optical pumping 
leads to a squeezed bath type coupling, and thus to symmetry breaking of 
decay on the Bloch sphere. With this system it should be possible to observe 
the effects predicted in the context of squeezed bath - atom interactions. 
The làser parameters allow one to choose properties of the squeezed bath 
interaction, such as the (effective) photon number expectation number N 

1 



and the squeezing phase <j>. We present results of a detailed analytical and 
numèrica! study. 
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I. INTRODUCTION 



The interaction of atòmic systems with squeezed light leads to interesting novel effects 
in atòmic dynamics |ï|. In particular, Gardiner || has shown that the atòmic Bloch vector 
of a two-level system coupled to a squeezed bath, which is characterized by a mean photon 
number N and squeezing parameter M, obeys the equations 

±(S X ) = - lx (S x ), (1) 

^(Sy) = -ly(Sy), 

±(S Z ) = - lz (S z ) - r 

The decay constants in this equation are given by 

7x = 7(iV+^-M), (2) 
7x = 7(iV+i + M), 
72 = 7 (2iV+l) 

with 7 is the spontaneous emission rate in free vacuum. We will refer to Eq. (JJ) as the 
Gardiner-Bloch Equation. According to Eqs. (0,0) the two quadrature components of the 
atòmic polarization (S x ) and (S y ) will decay with different rates ^ x and j y , respectively. In 
the limit of large photon number N ^> 1, and maximal squeezing, M 2 = N(N + 1), the decay 
of (S x ) will be suppressed according to 7 X — > j/(8N) (<^C 7 ) in comparison to spontaneous 
emission in free vacuum. At the same time the decay of (S y ) will be enhanced: 7 y — > 2^N 
(^> 7). As studied in numerous theoretical papers, these phase sensitive decay rates will 
also be visible in the spectrum of resonance fluorescence |3|,|4j], and the atòmic absorption 
spectrum M . A study in the context of the Jaynes-Cummings model has been done in M . 
From a physical point of view, the suppression of the decay below the free vacuum level is 
due to the reduced quantum fluctuations of one of the quadrature components of squeezed 
light. Similar effects have been investigated for the case of three-level systems For 
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experimental realization one has to include finite bandwidth effects. This has been done for 
example in P-p~3|j. 

On the experimental side, there have been only few experiments where the dynamics of 
àtoms in squeezed light has been studied in the laboratory. Most notable are the experi- 
ments by Kimble and coworkers who, for example, reported experimental observation of the 
linear intensity dependence 0JT^] of two-photon absorption rate in squeezed light jÏ6[ . The 
predictions of a suppressed and phase-sensitive decay of the atòmic polarization according 
to Eqs. (|Ï].H) have so far eluded observation, essentially due to the difficulty of embedding 
àtoms in a squeezed vacuum in a complete 4tt sòlid angle. If we denote by e the fraction 
of the sòlid angle filled by the squeezed vacuum modes then the polarization decay rates 
reduce to [01 



Ix = 7 
ly = 7 



e (JV + ±-M) + (l-e)± 
t (N+ l - + M) + {l-t) 1 - 



Clearly, if e differs significantly from 1 the influence of the squeezed vacuum will be reduced 
accordingly. One possible solution to achieve a large effective e close to 1 is to consider 
systems which are effectively one-dimensional due, for example, to a strong cavity-atom 



interactions as proposed (see, e.g. Refs. [[Ï7j and fïï 

In this paper we will show how to "engineer" a squeezed-bath type interaction leading 
to a Gardiner-Bloch Equation ([!]) by studying the dynamics of a driven multilevel atom 
coupled to " normal " vacuum. In the simplest version we will consider a four-level system: 
for example, an angular momentum J g = 1/2 — > J e = 1/2 transition with two (degenerate) 
ground (|<7 m =±i/2)) an d excited states (|e m= ±i/ 2 )), as illustrated in Fig. |j. If this atòmic 
transition is driven by a + and a_ polarized làser light, the spontaneously emitted linearly 
polarized 7r-photons emitted in the transition |e + i/ 2 ) —> [(7+1/2) and |e_i/ 2 ) —> IP-1/2) will 
interfere since they are indistinguishable. For weak driving fields far below saturation we 
can adiabatically eliminate the excited states. The dynamics of the two ground states 
| í7 + i/ 2 ) , |<7-i/2) then obeys a Master equation with damping terms due to optical pumping 



processes between the two ground states. This master equation has a structure analogous 
to the coupling of a two-level àtoms to a squeezed bath. ([I]) Interfering processes to mimic 
squeezed state detection statistics have been used previously by Wilkens, Lewenstein and 



Grochmalicki. |i9| , |20| It is used in the context of atòmic spin measurements in |2Ï[ . Reservoir 
engineering to influence resonance fluorescence by changing the density of states in a cavity 
has been investigated in [p2|. In the context of ion motion reservoir engineering has been 
shown in [ 23(1 . 

The paper is organized as follows: In the section two we will summarize the properties 
of the master equation for a two-level atom coupled to a squeezed bath for reference and 
comparison in later sections. In section three we discuss the reduction of the multilevel 
master equation to an effective two-level master equation with squeezed-bath-type couplings 
of the form (JJ). Sec. five presents numerical results for resonance fluorescence and absorption 
spectra in four level systems, and compares with the corresponding squeezed bath results. 
Non-ideal effects, leading to thermal (as opposed to squeezed) bath coupling, as given for 
example by the cross decay terms in spontaneous emission in Fig. [I], will be investigated in 
Sec. 6; in addition we will suggest mechanisms how to suppress these unwanted effects. 



II. TWO LEVEL ATOM IN A SQUEEZED VACUUM: A SUMMARY 

In this section we will briefly review the bàsic effects and properties of a two-level system 
coupled to a squeezed bath. In particular we discuss solutions of the Gardiner-Bloch Equa- 
tions JJ, and the main features of the spectrum of resonance fluorescence and the atòmic 
absorption spectrum. We summarize these results for reference in the following sections. 



A. The Gardiner-Bloch Equation 

We consider a two-level atom with ground state \g) and excited state |e). The coupling 
of the atom to a squeezed vacuum is described by the interaction Hamiltonian 
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H mt = art + a t r . ( 3 ) 

with atòmic lowering operator is a = \g){e\ and a bath operator T, which given in terms of 
the coupling constants A and the photon annihilation operator a^A as 

T(t)=J2^,xa k , x e- iuJkt . (4) 

k,X 

A broadband squeezed bath centered around the atòmic transition frequency uja is charac- 
terized by the correlation functions 

(^t(í)^(í')> = 7iV5(í-í , ), ( 5 ) 

(T(t)T^t'))= 1 (N + l)5(t-t'), 
(r(í)r(í')) = ~iMe-' l,t, e- 2ïu)At 5{t - t'), 
(r t (í)r t (í / )) = ^Mé*e 2iuJAt 5(t - t'), 

with effective photon number N and the (real) squeezing parameter M. These parameters 
are restricted by the inequality M 2 < N(N + 1), where the equal sign holds for maximal 
squeezing. The squeezing phase is denoted by <fi. Here "broadband squeezed bath" refers 
to the assumption that the squeezing bandwidth is larger than the other frequency scales in 
the problem (excluding the optical frequency), such as the spontaneous decay rate 7. 

Knowledge of the correlation functions (|5[) allows one to derive the master equation for 
the atòmic dynamics in the Born-Markov approximation. In a rotating frame the master 
equation is m 



d 

dt r 2 



p = I 7 (iV + 1) (2apa ] - aVp- ptrV) (6) 
+ \lN (2c 



_ } . ^-Za^pa — acr' p — paa^ 

For non-perfect squeezing we define mean photon numbers N\ and N2 through Ni(N\ + 1) 
M 2 and N = Ni + N 2 , which allows us to rewrite (|6]) in the form 
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= ~ 7 (2SpSt - SVS - pSts) (7) 



+-7ÍV2 {2a pa^ — a^crp — pa^a 
+ -7ÍV2 (2o^ po — <7<7^ p — poo^ 



dt r 2 

2 
1 

2 

with 

E = 0VÏ + la + e^^ïvï^. 

The first line in Eq. (^) are damping terms for ideal squeezing, while the second and third 
line in correspond to a thermal reservoir (background) . 

In a wave function simulation of the master equation we would intèrpret E, and a, cr^ as 
quantum jump operators for the various damping terms in (0). In this quantum jump picture 
the master equation is represented by an ensemble of quantum trajectòries of pure system 
wave functions, where the evolution of states is described by an effective, non-hermitian, 
Hamiltonian H e s , interrupted by quantum jumps where the wave function undergoes jumps 
according to the action of quantum jump operators on the wave function. For a discussion 



of a quantum jump picture of two-level àtoms coupled to a squeezed bath we refer to p4 

With the notation S x = ^ + cr^>, S y = Mcr^ — en /i}, S z = (P e — P g ) for the Blochvec- 
tor, the Master equation (||) with the phase choice (p = is equivalent to equations (D) given 
in the introduction. The corresponding equation for a driven system reads 

^(S x ) = -i{N+ l - -Mcoa<f>)(S x ) + 7Msin0(5,}, (8) 



—À S v) = + - + M cos 0) + 7Msin0(S x > + Ü D (S Z ), 



±(S.) = -j(2N + l)(S z ) - n D (s y ) - 7. 

where <fi = <p g — 2(p d with cpu as phase of the driving field with Rabi frequency and (pg 
as squeezing phase in the same reference frame as (pr>. The Bloch equation is given here 
for the choice (po — 0. A consequence of the broken symmetry in the polarization decay 
is that the steady state of the driven two-level system becomes dependent on the relative 
phase between the squeezing and the driving field. 



B. Steady State 



The steady state of the Bloch equation (|j) for the driven two-level system becomes 
phase dependent if the vacuum is squeezed (M ^ 0). This means that propagation effects 
like absorption and dispersion become dependent on the the relative phase between the 
driving field and the squeezing phase. Therefore the steady state might be interesting to be 
used as a kind of optical switching where propagating beams are controlled by the squeezing 
phase and the driving làser. For normal vacuum the steady state is (up to a global phase) 
independent of the phase of the driving làser. 

C. Spectrum of Resonance Fluorescence 

The spectrum of resonance fluorescence spectrum of a strongly driven two-level atom is 
determined by the Fourier-transform S(u) of the stationary two-time correlation function 
of the atòmic dipole: 

S 2 {u) = FT(J{0)<t{t)). (9) 

Eq. (R|) implies that the free field corresponding to the particular mode for which the 
spectrum is observed is not squeezed. 

The spectrum of resonance fluorescence for a squeezed bath has been calculated in Ref. 
||. According to the Quantum Regression Theorem the evolution of the correlation function 
is governed by the Bloch equations. Therefore the spectral lines of the resonance fluores- 
cence spectrum are determined by the eigenvalues of the coefficient matrix in (§). For strong 
driving ^> 7 X ·,7 Ï/ one finds a Mollow triplet with phase dependent line-widths and inten- 
sities. (See table |[) For strong squeezing (N — > oo) the center line can become arbitrarily 
close to zero line-width for = 0. 
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D. Atòmic Absorption Spectrum 



The absorption spectrum of a weak probe field in presence of a strong driving field is 
related to the Fourier Transform of the two-time correlation function 

W{u) = FT{[<t(t),<t^0)]). (10) 

According to Ref. || W(u) shows strong dependence on the value of (fi. For (fi = it shows a 
strong peak at the center and dispersion-like sidebands. The width of the central absorption 
peak is again dependent on N and narrows to zero as N grows since its linewidth is the 
same as in the resonance fluorescence. Between the center line and the sidebands is a region 
with gain. For (fi = ir these detuning regions with gain are stronger and one finds even for 
the center line stimulated emission. The center line is broadened as well in accordance with 
the center linewidth of the resonance fluorescence. 

III. DERIVATION OF A SQUEEZED BATH MASTER EQUATION FOR AN 

EFFECTIVE TWO-LEVEL SYSTEM 

The signature of the squeezed bath coupling in the master equation (|?p is the appearance 
of a quantum jump operator of the form y/Ni + la+e^y/T^a^ . We intèrpret the structure of 
this operator as arising from interference in the atòmic transition from the lower to the upper 
atòmic level (ert), an( j th e transition from the upper to the lower state (er). An analogous 
interference will oceur in optical pumping processes between the two ground states in a 
four-level system in the emission of 7r-polarized photons. 

Let us consider the four-level atom according to Fig. |I[ The two ground states \g m =±i/2) 
correspond to the two levels \g), |e) of the previous section. The two upper levels |e m= ±i/ 2 ) 
will be used as auxiliary levels. The upper and lower levels are connected by spontaneous 
decay due to the interaction with the normal vacuum reservoir. The two decay channels 
are assumed to give the same decay rate T and to give rise to emitted photons of the same 
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polarization and frequency, so that the two spontaneous emission processes are indistin- 
guishable. We add two weak resonant làser fields to connect col·ierently the levels lflL-1/2) to 
le+1/2) with a Rabi frequency e_í2 (<C T) and | (7+1/2) to |e_i/ 2 ) with a Rabi frequency e+í2 
(<C T). (e + + = 1) Thus, the two ground states are connected by optical pumping. A first 
process starts at level | (7-1/2), transfers via the weak làser field with Rabi frequency e_f2 to 
level le+1/2) from where it decays to level | (7+1/2)- The second process, will connect | (7+1/2) to 
|<7-i/2)j due to absorption from the other làser with Rabi frequency e + Q (<C T). As a result 
we end up with a master equation for the effective two-level system | (7+1/2), \g~ 1/2) where the 
damping terms have the structure analogous to a squeezed bath coupling. The cross decay 
terms corresponding emission of a cr+ or cr_ photon will not interfere, and thus give rise a 
phase-insensitive background (deviation from the ideal squeezed-bath type couplings). 

A. Adiabatic Elimination in the Four-Level System Master Equation 

We will now derive the master equation for the effective two-level system starting with 
the master equation for the four-level system, including the cross decay (Fig. |ï|) It is given 
by 

^p = -z (H ca p-pHl s ) (11) 
+gfT {at + (x 2 ) p (a\ + 4) 

with H e g an effective Hamiltonian 

H cS = -i^P e + e_| (VI + o-_) + e+| + é^a + ) (12) 

where <f>L is the relative phase between the two làsers. The atòmic lowering operators are 

fi = 1 0-1/2) (e_i/ 2 |, o"2 = |#+i/2)(e+i/ 2 |, 
°"+ = l5 , +i/2)(e-i/ 2 |, <r_ = | g-y 2 ) (e+i/ 2 | • 
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Furthermore, we define a projection operator for the upper states as 

P e = |e_i/ 2 )(e_i/ 2 | + |e + i/ 2 )(e + i/ 2 |. 

The Clebsch-Gordan coeficients for the coupling of the upper to the lower levels are denoted 
by gi and g c , respectively. (Compare Fig. []].) 

In a parameter regime satisfying e±íl C T we can adiabatically eliminate the upper 
levels |e_i/ 2 ) and |e +1 / 2 ) • For this we introduce the operator 



D = e_a_ + e+e^a. 



-i 



and use the projections onto the upper states, P e , and the lower states, P g , to represent 
the master equation in terms of the evolution of the upper states, the lower states and the 
coherence between upper and lower states. The corresponding equations are 

^P e pP e = ^p ee = ~Tp ee - i^ (p ] p ge - PegD) (13) 
^PgpPe = ^Pge = ~Pge + Í~ {pggD ~ Dp ee ) (14) 
J t P 9P P H = J t P99 = ( D Peg - PgeD ] ) (15) 

+gfT (cri + cr 2 ) p ee (a\ + a]) + glYa^p ee aL + g 2 c Ta + p ee a\_ . 



In the adiabatic approximation, p ee follows adiabatically the changes in the coherences de- 
scribed by p eg which in turn follow adiabatically the dynamics of the ground states p gg . In 
this approximation we find for the density matrix restricted to the lower levels, p gg , a master 
equation which involves two Lindblad terms, 

^ t p 99 = \y 9 " ( 2èpèt " t]tp - pt]t ) (16) 

+ l 9 *T {2azpcTz ~ 2p) ' (17) 

with the jump operator S given in terms of the Raman operator o = |<7_i/ 2 )(<7 + i/ 2 | as 

S = (e + a + é*e_J) . (18) 

The remaining parts of the four level density matrix are related to p gg by 
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> ee = %Dt Paa D (19) 



and 

Peg = i^Dt Pgg . (20) 
We make the replacements of table |T| to write this master equation as 

= -i (H cS p - pHl s ) + (21) 

fe 



2 l'N 



with the effective Hamiltonian 

H,n = -i\ [gl ( f + i) + 3? (ml + P+)) • (22) 

The additional jump operator damping term involving 

v z = P + - P- 

is characteristic for a process which destroys the decoherence between the two ground state 



levels. [£(J Now we also consider a làser connecting the Raman transition | g-1/2 ) — * | (7+1/2)- 
We will assume that these làsers are far from any atòmic resonance so that spontaneous 
emission can be neglected. (Alternatively, one can consider an AC magnètic field.) The the 
Bloch equation is given by 

±(S X ) = -j{N+ X - -gfM cos <p L )(S x ) + 19 ?M sin <p L (S y ), (23) 
f t (S y ) = -j{N + i + g?M cos <j) L )(S y ) +jgfMsm<t> L (S x ) +iï D (S z ), 
±(S.) = -igf(2N + 1)(S Z ) - íi D (S y ) - gh 

where is the Rabi frequeney of the Raman transition with the phase chosen as <$r = 0. 
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B. Discussion 



For gi = 1 (no cross decay) Eq. ( p3[) is precisely the Bloch equation (|8]). We can bring the 
master equation (|i6|) into the form (||) or the Bloch form ([]]) by making the identifications 
of table O. The two ground levels thus show a dynamics which is analogous to that of a 
two-level atom in a squeezed bath. Effectively we "engineer" with help of the upper levels 
and the làser fields a a reservoir for the two ground levels such that it looks like the coupling 
to a squeezed bath. (See figure This point of view can be brought out more clearly by 
writing an interaction Hamiltonian for the lower states in the adiabatic elimination as 

#int = Vggï f + aj g f (24) 

with a gg = | (71/2) (fl 1 - 1/2 1 the transition operator between the two ground states, and the 
bath operator V given in terms of the coupling constants Kk,\ and the photon annihilation 
operator a& \ as 



E 

A-.A 



1- l i« fciA e_aí ^e" 1 ^* - < ^e^a^e^e^ 



(25) 



The correlation function of this bath operators are the same as the correlations functions of 
the operator T given in (||) of the squeezed bath if one uses the identification of table [TT] and 
takes into account that for the degenerate ground-states the frequency corresponding to u>a 
vanishes. 

In summary, the main elements of the four level system is that the direct transition 
between the two levels of the two-level scheme is replaced by a pumping process involving 
other atòmic levels, and that the driving between the two ground states is replaced by a 
Raman coupling. 



IV. DISCUSSION AND RESULTS: THE IDEAL MODEL 



We will now present some examples of the ideal (gi = 1) behavior of the four level 
system to illustrate that it indeed mimics a two-level atom coupled to a squeezed bath. The 
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following figures were produced using the full master equation for the four level system of 
Fig. |ï| Analytic formulas presented, however, use the adiabatic eliminated equations for 
the ground states to allow an easier comparison to the two-level system. For the parameters 
used in the figures, the adiabatically eliminated and the exact calculations are in excellent 
agreement. 

A. Steady State Solution 

An example of exact correspondence is the steady state solution for the driven system 
as investigated by Carmichael, Lane and Walls. We drive the system by a resonant 
Raman transition with an effective Rabi frequency Qp. A typical steady state is shown as 
sòlid line in figure [| The dependence on the phase <fi allows light fields, causing the two 
internal transitions and the driving, to interact depending on the relative phase between 
them. This leads to effects in the propagation of those fields through a cloud of àtoms with 
the effective four level system as shown above. The results for the four-level system show 
the 0-dependence in agreement with a corresponding two-level system. 

B. Absorption Spectrum 

As our first example we study the absorption spectrum of the four level system. The 
idea is, in analogy of the investigation concerning the two-level system, to drive the effective 
two-level system strongly coupling the two ground states in a Raman transition which is 
tuned on resonance. Then the weak field absorption spectrum is measured, where the probe 
field is represented by another stimulated Raman transition. The stationary absorption 
spectrum is then given by the Fourier transform of a two-time correlation function 

W(u) ~ FT(a gg (r)al g (0) - al g (0)a gg (r)) (26) 

which can be calculated using the Quantum Regression Theorem. 
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The absorption spectrum of the four-level system (see figure |] (a) ) is identical to that 
of the two-level system in a squeezed bath and we can see the phase dependence and the 
sharp peak at the center for = (sòlid line). This is readily understood since the process 
involved in measuring the absorption spectrum makes use only of the system level dynamics 
as described by the master equation (21). This stands in contrast to the situation of the 
resonance fluorescence which differs slightly from that of the two-level scheme, as will be 
shown in the next subsection. 



C. Resonance Fluorescence 

For resonance fluorescence we drive the system again by a stimulated Raman transition 
and detect the fluorescence coming from the transitions \e±i/2) — > \g±i/2)- This leads, in the 
adiabatic elimination, to the spectrum given in terms of ground state correlations as 

S(u>) = FT ((e+MjO) + e.a gg (0)) (e + e^a gg (r) + e-Jjrj)) . (27) 

In contrast, the fluorescence spectrum calculated by is based on the Fourier Transform 
of the correlation function S 2 (uj) = FT ^ert(0)cr(r)^ (see Sec. 2), which is based on the 
assumption that there is no interference in the detector between the source field and the 
squeezed vacuum modes. In our case the fluorescence spectrum is due to a kind of atòmic 
quadrature correlations (similar to those appearing in the squeezing spectrum ||271| ) and 
there are three additional terms with respect to the two-level scheme if we use the four level 
scheme. However, this does not change the position and width of the lines in the Mollow 
triplet. Therefore the line narrowing and the phase dependence of the line intensities can be 
observed in the four-level system. (See figure |^ (a).) The resonance fluorescence of the four- 
level scheme differs from that of the two-level which is not surprising since its origin of the 
can be explained only in the full four-level scheme and not in the effective two-level scheme 
of the groundstates. In this point the resonance fluorescence differs from the absorption 
spectrum and the steady state which can be explained just in terms of the Master equation 
for the two ground states. 

15 



We continue to speak of the Mollow triplet for the four-level system in this parameter 
regime since only these three lines have non-negligible intensity. As one leaves the regime of 
validity of the adiabatic elimination, where the effective two-level system and the four-level 
system are in perfect agreement, one finds, of course, a richer line structure. (See subsection 
VTJ .) 



V. RESULTS AND DISCUSSION: NON-IDEAL EFFECTS 

In this section we will address separately two problems connected to the validity of the 
effective two-level system master equation with squeezed-bath type couplings. One is the 
question of the influence of dephasing terms in the master equation fl2lD due to cross decay. 
The other is the question of behavior of the system at the onset of saturation. 

A. Effects of Cross Decay 

For a four-level system J g = | — ► J e = \ the Clebsch-Gordan coemcients are gf = | 
and gl — |, so that most of the spontaneous processes go into the (unwanted) dephasing 
part. This ratio improves by using Zeeman sub-levels of higher angular momentum as upper 
levels where with help of the AC Stark effect only the levels with m = ±| take part in the 
dynamics while the other levels are shifted off resonance by a làser field. For J e = | we then 
find the slightly more favorable numbers of gf = | and g\ — g. 

The effect of a non-vanishing g c on the spectra and the steady state is quite different. 
The absorption spectrum and the resonance fluorescence spectrum show even for g% = | a 
strong dependence on 0. (See figures [| and ||.) 

The line-width of the resonance fluorescence interpolates linearly as a function of g\ 
between the vàlues vàlid for the interaction with a squeezed vacuum and those for the 
normal vacuum as shown in table |iTT| . The steady state variation with <f), especially that 
of the component S x , is more sensitive to the presence of cross decay. To utilize the <fi 
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dependence of the steady state one has to find means to suppress this cross decay to an 
extent as large as possible (see figure |3|.) 



B. Suppressing the Effects of Cross Decay 

To suppress the dephasing in favor of the squeezed-bath-type effects one can use de- 
structive interference of the cross decay in a configuration with more atòmic levels, similar 



to the one proposed in Pq| . This scheme employs upper levels Zeeman levels of angular 
momentum J e = ~ and those of angular momentum J a = |. (See figure |5|.) The weak làser 
fields of frequency u L are now detuned between the upper levels with J e = ~ and those with 
Ja = \ with detuning A e = ujl — uj eg and A a = — (u>l — uj ag ), respectively. The spontaneous 
decay rate and the Clebsch-Gordan coeflicients have subscript indicating the level they are 
referring to. We can calculate the rate for processes starting and ending in the same ground 
level. This rate depends on the detuning and the dipole elements d e and d a , referring to the 
transition between the selected ground state and the two intermediate upper levels taking 
part in the transition respectively. It is given by 

1 12c 3 e |£| 



T hu) 3 



L 



(») 2 r ( e ) 2 r 

í/c L a y c L e 



A _ iEa A _ íík 

a 2 e 2 



(28) 



where we introduced the field strength \E\ corresponding to the Rabi freqüències e + /_í7 

2 

... "r a 

the cross decay by an order of magnitude in jr^-. A demonstration of this effect in a similar 

^e/a 



respectively. We always can choose the detuning such that ^ = 9c (a) l e . This will reduce 



scheme has been given by Xia, Ye and Zhu |29|]. Any remaining cross decay can be described, 
within the validity of the adiabatic elimination, by effective Clebsch-Gordan coeflicients and 
an effective decay rate. [EÏÏJ 



C. Validity of the Adiabatic Elimination 

Finally, we investigate the behavior of the four-level system if we increase the Rabi 
freqüències of the infernal transitions. To explore the range of validity of the adiabatic 
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elimination and the effects of onset of saturation, we compare line-width and position of the 
Mollow triplet for the adiabatically eliminated situation and the full four level system. In 
the latter one we characterize the Mollow triplet as the three narrowest lines in the resonance 
fluorescence spectrum. We can find their positions and line-width as eigenvalues of the 15 
by 15 matrix appearing in the Bloch equation of the four level system. The effect of an 
increased value of Q/T for the parameters used to plot the resonance fluorescence in figure 
H] is shown in Fig. [7[ Up to Q/T rs 0.2 the line width and positions of the two-level Mollow 
triplet follow the curve predicted by the adiabatically eliminated theory. For higher vàlues 
of Q/T the same position is still predicted by both equations but all three lines become 
narrower in the full system than expected by the reduced equations. Other lines of the 
full four level system prové to be clearly distinct in their line-width from the ones of the 
Mollow triplet. We have checked that the whole appearance of the resonance fluorescence 
(that is including the line amplitudes) shows an increasing deviation from the adiabatically 
eliminated equations as Q/T increases beyond the value Q/T = 0.2 just as expected from 
the behavior of the eigenvalues. 



VI. CONCLUSIONS 

The central result of this paper is that there is a way to engineer an environment for 
a two-level system. The bàsic idea is the use the regime of adiabatic elimination to link 
up coherent transitions from system levels to additional atòmic levels with spontaneous 
transitions from these additional levels back so system levels. This leads to jump opera- 
tors in the master equation which can be designed to purpose. Especially interesting jump 
operators can be achieved if one uses, as in the model presented here, that there are indis- 
tinguishable spontaneous processes leading to linear superpositions of jump operators. With 
this engineering one converts the trivial reservoir of normal vacuum modes into the more 
sophisticated one of an effective squeezed vacuum filling the whole sòlid angle of 4ir. 

We have illustrated this procedure for the example of a squeezed bath-type coupling. 
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Within this example we demonstrated this reservoir can be engineered by the use of Zee- 
man sub-levels. The predicted effects of the squeezed bath coupling on the steady state, 
the absorptions spectrum and, with a slight change, the resonance fluorescence have been 
recovered in this model. We have shown that one can select spontaneous decay channels 
in the multi-level scheme and suppress others. Therefore one can enhance the effect of en- 
gineered jump operators in the final master equation and suppress other processes which 
would partially destroy the desired effects. In a numerical study we have shown the effects of 
a non-ideal realization resulting in decoherence processes which are collision like. This study 
shows that the predicted effects is robust enough to allow their observation. The regime of 
validity of the adiabatic elimination of the upper levels in our model has been explored. 

Our work presents an approach which allows the observation of theoretically predicted 
effects within reach of todays experimental techniques. It encourages the investigation of 
other reservoir couplings which could be engineered with the ideas presented in this paper. 



This work has been supported by the by the Austrian Science Foundation and by the 
EU as part of the TMR-Network "Microlaser and Cavity QED" ERB-4061-P1-95-1021. 



The value of the detuning which gives destructive interference can be found using second 
order perturbation theory. We consider the subsystem as shown in figure || The Hamiltonian 
describing the subsystem is given by 
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APPENDIX: TRANSITION RATE FOR CROSS DECAY 



H = H A + H B + H AB + H AL 



(Al) 



with the system Hamiltonians of the atom, the bath and the làser mode given by 



H A = k(j e P t 



e- 1 e+i/2 



a A 0+1/2 



(A2) 



19 



fcA 



and the interaction Hamiltonians atom-bath and atom-laser given by 
H A b = [ hK k\ a k\\e + i/2)(g-i/2\ + hn^a k x\2a)(g_ 1/2 \ + h.c. 

fcA 

+ Y. [^4A a fcA|e-i/2)(5'+i/2| + l a k\ I la) (0+1/2 1 + h.c. 

fcA 

H A L = hQ( e )\e +1 / 2 )(g-i/2\ + hQ( a )\a +1/2 )(g-i/2\ + h.c. 



+ 



(A3) 



(A4) 

(A5) 
(A6) 



Here the K kX , fcjj." are coupling constants defined in terms of dipole vectors d( e / a ), unit vectors 
for the elèctric field e^x and a quantization volume V by 



K 



(e/a) 
fcA 



í^fcA 



(e/a), 



(A7) 



V 2^e y 

The Rabi freqüències due to the interaction with the classically treated làser fields with 
elèctric field of amplitude \E\ and unit vector ei are 

2\E\ 



n 



(e/a) 



{&L-d(e/a)) 



(A8) 



In second order perturbation theory the transition rate for processes starting and ending in 
\g-i/2) under spontaneous emission of a circular polarized photon is given by 



1 A\E\ 2 uj\ 



\d a 



\d,, 



T 37TC 3 h 3 € 

12ch \E\ 2 
hu 3 L 



A a A 

9',"" 



r. L a y r . L e 



(A9) 
(A10) 



Inclusion of the main terms of higher orders is by the resolvent method pTJ yields the result 



1 _ 12c 3 e |£| 

T 



o (a)2 r 

i/c 1 o 



A — i^-a- A — 

a 2 e 2 



(All) 



This leads finally to the expression (p8|) . The same result can be obtained using adiabatic 
elimination in the master equation for the full system as shown in figure || 
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TABLES 





= 


(j) = 7T 


central peak 


Ix = 


(7y + 7*)/2 = 




7 (N + \ - M) 


7 (iV + \ + Af) 


sidebands 


% = 


(7x+7.)/2 = 




l (6N + 3 + 2M) 


|(6A + 3 - 2M) 



TABLE I. Line-width in the regime of strong driving. The notations Jxtlyjlz refers to the 
notation of equations (Q). 



two-level atom in squeezed vacuum 



four level system 



spontaneous emission rate in 
normal vacuum 
photon number expectation 
value of squeezed vacuum 
squeezing parameter 
phase of squeezing 




(4 - £) ¥ 



72— T 



72 ~T 



=1<PL 



TABLE II. The parameters of the two-level system in a squeezed vacuum expressed as a 
function of the parameters of the mimicking 4 level system 





<p = 


4> = vr 


central peak 
sidebands 


j(N + ±- gi 2 ^/N(N + lj) 

l (2N + 1 + 2 9l 2 {2N + 1 + y/N(N + 1))) 


y(N+^ + gi 2 y/N{N + l)j 
l(2N + l + 2 9l 2 (2N + l- 





TABLE III. Line-width in the regime of strong driving in dependence of effective photon 
number N and Clebsch Gordan Coemcient gf = 1 — g 2 . 
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FIGURES 

FIG. 1. Realization of a level scheme involving two ground states and two upper states which 
are Zeeman sublevels. The atòmic levels are coupled by two right and left circular polarized làser 
fields with Rabi freqüències e+fi and e_í7. These làser fields are weak in the sense that T ^> e±íl. 
We will consider this scheme in two situations: (A) In the ideal situation there is only spontaneous 
decay along channels with Anij = 0. These decays interfere and give rise to squeezed bath like 
effects. (B) In a realistic model each upper level decays with a total decay rate T with a branching 
probability for the decay channels determined by the Clebsch-Gordan coefhcients g c and g\. The 
processes along the cross lines of decay do not interference since the give rise to right and left 
circular polarized photons respectively. The effect of these cross-decay channels will give rise to 



collision-like effects of the reservoir. In subsection V B we show how to suppress the cross-decay. 



FIG 2. Schematic view of the relation system to reservoir for the two-level and the four level 
situation. 

FIG 3. Steady state for N = 2.1, J = 1.9 x 1(T 5 , = 5.1 and for the ideal case of 
gi = 1 (sòlid line) and three other vàlues with gi = 0.99 (dashed line), gi = 0.95 (dash-dotted line), 
gi = 0.9 (dotted line). 

FIG 4. Absorption spectrum (in arbitrary units) for N = 1, ^ = | x 10 -4 and ^p- = 7.1 and 

/ „\l/2 /iX 1 / 2 

for (a) the ideal case (gi = 1) and for three vàlues of g\ = {0.9 , ( 3 ) > ( 3 ) } (sub-plots (b) to 
(d)). The sòlid (dotted) lines show the spectra for = (<fr = n). 



FIG. 5. Resonance fluorescence in arbitrary units for N = 0.2, ^ = 7.1 x 10 5 and — = 7.1 

1/2 



/„\ i/z / 1 \ í/Z 

and for (a) the ideal case (gi = 1) and for three vàlues of g\ = {0.9 , ( 3) > ( è) } (sub-plots 
(a) to (d)). The sòlid (dotted) line shows the spectrum for <fi = (<f> = 7r). 

FIG 6. Extended level scheme which is used to suppress spontaneous decay between levels 
with mj = — \ and mj = +\ (cross decay) due to destructive interference. 
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FIG. 7. Study of line positions and width for the full four level system (sòlid and dash-dotted 
lines) in comparison with the results of adiabatically eliminated equations (dotted lines). The 
parameter vàlues are N = 0.2, 2a = 7.1 and g x = 1. Compared are the line-width of the center 
line (top left), the line-width of the sideband lines (bottom left), the positions of the three lines 
(top right). To demonstrate that the Mollow triplet is clearly distinguishable from other lines we 
show the line-width of the line which is the narrowest after the triplet. (Bottom right.) From this 
comparison it is clear that for Íï/T < 0.2 the approximation made in the adiabatic elimination of 
the upper states is vàlid. 

FIG. 8. Subsystem relevant to investigate the destructive interference of spontaneous emission 
along the pumping transitions. The total decay rates of the upper levels are T e and r a . The 
Clebsch-Gordan coefficients for the spontaneous decay emitting circular (linear) polarized photons 
are denoted by and gf^ (g\ e ^ and g^). The làser is described by the elèctric field amplitude 
\E\ and shown as sòlid arrow. 
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